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' We analyze the scattering sector of the Hamiltonians for both gapless and gapped graphene in 

(~| , the presence of a charge impurity using the 2D Dirac equation, which is applicable in the long 

^ • wavelength limit. We show that for certain range of the system parameters, the combined effect of 

^ ' the short range interactions due to the charge impurity can be modelled using a single real parameter 

appearing in the boundary conditions. The phase shifts and the scattering matrix depend explicitly 
on this parameter. We argue that this parameter for graphene can be fixed empirically, through 
measurements of observables that depend on the scattering data. 
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INTRODUCTION 



^ ' The experimental fabrication of monolayer graphene [H, 0] , a two dimensional system with a hexagonal honeycomb 
structure, has opened up an area of physics which is of wide theoretical and experimental interest. The low energy 
fN| . dynamics of gapless graphene is described by a 2D massless Dirac equation [sj. The effect of charge impurities in 
' gapless graphene has been studied extensively, where bound states are predicted when the Coulomb charge exceeds 
CO a certain critical value 0, 0,H0]- Below this critical value, the bound states do not form in massless graphene due 
to the Klein paradox. The generalization of this model where a Dirac mass term is included in the Hamiltonian 
. is known to describe gapped graphene. The effect of charge impurities in gapped graphene have also been studied 
J> ! [li S 0j [m I where bound states can form for arbitrarily small value of the Coulomb charge. 

' Apart from providing an axially symmetric Coulomb potential, the charge impurities may also induce other short- 
5^ • range or singular interactions, such as a delta function type potential. The detailed knowledge of such short range 
. 5^ interactions in any given sample would be very hard to ascertain empirically. Moreover, the 2D Dirac description for 
graphene is valid only in the low energy (long distance) limit and it would not be practical to include short range 
interactions directly in such a Hamiltonian. One possible way to model the combined effect of such short range 
interactions on the long wavelength dynamics is through the choice of appropriate boundary conditions. 

Recently we have analyzed the effects of generalized boundary conditions, that may arise in the presence of a 
charge impurity, for the bound state sectors of gapless 0] and gapped graphene. In these papers, the choice of 
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appropriate boundary conditions due to the presence of charge impurities in graphene was guided by the principle of 
self-adjointness of the corresponding Hamiltonian. This is a well known approach that provides a reliable description 
of bound states in systems with delta function type potentials 13 1 as well as in fermionic 0, 0, [l^ E3i llS] and 



anyonic systems 19|, |20[ and in molecular physics [21]. It has already been used to study certain topological defects in 
graphene 2^ 27 1. Apart from bound states, this method is also useful to study the effects of boundary conditions on 
scattering states, phase shifts and S matrices in systems with short range interactions [22. 23, 24, 2^. Application of 



introduced a new parameter 



this approach in the analysis of bound states in graphene with charge impurities 0, [ 
which labels the boundary conditions. This self-adjoint extension parameter has to be fixed empirically, possibly 
through the scanning tunneling microscopy (STM) measurements of the local density of states (LDOS). 

In this paper, we analyze the scattering states of gapless and gapped graphene within the framework of an effective 
low energy description in the presence of a charge impurity. Here again we model the combined effect of the short 
range interactions due to the charge impurity through the choice of boundary conditions, which is guided by the 
principle of self-adjointness. In the scattering sector, the phase shifts and the S matrices depend explicitly on the 
boundary conditions, through the choice of the self-adjoint extension parameter. The knowledge of the phase shifts 
can be used to calculate physical properties of the system, such as resistivity 28|. In addition, the effect of this new 
class of boundary conditions can also appear in the density of states and Friedel sum rule. In our formulation, all 
these observable effects would depend on the self-adjoint extension parameter, which can be used to determine it 
empirically. For consistency, we also demonstrate that the bound state energies calculated in lj| can be recovered 
from the poles of the scattering matrix as well. 

This paper is organized as follows. In Section 2 we discuss the effect of boundary conditions on the scattering sector 
of the massless graphene with a charge impurity. In Section 3 a similar analysis is carried out for the gapped graphene 
with a charge impurity. We conclude this paper in Section 4 with a discussion of the experimentally observable effects 
of our analysis. 



2. SCATTERING IN GAPLESS GRAPHENE WITH CHARGE IMPURITY 

In this Section we shall analyze the scattering states, phase shifts and the S matrix for gapless graphene in the 
presence of a charge impurity. We begin by a brief review of the gapless graphene system when a charge impurity is 
present. Next we discuss the allowed class of boundary conditions for this system and end by doing the analysis for 
the scattering sector using the new boundary conditions. We show that the resulting S matrix and the phase shifts 
depend explicitly on a parameter which labels the boundary conditions. We argue that this parameter cannot be 
determined from theoretical analysis but must be fixed empirically. 

2.1. Gapless graphene with a charge impurity 

The 2D Dirac equation describing a Coulomb type charge impurity in a gapless graphene monolayer can be written 

as 

H^^E^, (1) 

with the Hamiltonian H being given by 

H = -thvFiaA + <J2dy) - ^ ^ hvp -'^-^ ) - 7> (2) 

where r is the radial coordinate on the two dimensional x — y plane. In above, vp is the graphene Fermi velocity and 

y 2 

a = where Z denotes the impurity valence, e is the unit charge and k is the effective dielectric constant. One 
can separate Eq. ([T]) by assuming the solution of the form 

/ F{r) 1 e*0-^)'> \ 
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where (f) denotes the corresponding polar angle and j is the half integral azimuthal quantum number. Consider the 
ansatz 

F(r) = e''=V^-5(u(r) + w(r)), (4) 

G(r) = e'''W~^u{r) - v{r)), (5) 
where v = J — (P- , B = -r^ = f ° , k = —-iM—. This leads to the radial Dirac operator 

and the set of coupled equations 

du 

r— + {v + iP + 2ikr)u- jv = 0, (7) 
dv 

r— + (i/ - il3)v - ju = 0. (8) 

After introducing a variable z = —2ikr, the last two equations can be combined into a single equation for the 
function v, given by 



z^ + {l + 2iy- z)^-{iy-il3)v = 0. (9) 
az^ az 



This shows that v is described by a confluent hypergeometric function 29| and a similar conclusion holds for the 
function u as well. 

In what follows, we shall set h = vp = 1. It may be noted that since the minimum magnitude of the azimuthal 
quantum number is j = ^, the corresponding critical value of /3 is given by f3crit = \- In this paper, as far as the 
massless graphene case is considered, we shall be interested only in the subcritical region. In the subcritical region, 
the parameter v is a. positive real number with v = Q denoting the critical value. 



2.2. Boundary conditions for gapless graphene with a charge impurity 

Equation ^ defines the radial Dirac operator H^. for a gapless graphene with a charge impurity located at r — Q. 
As discussed before, the impurity can give rise to additional short range interactions which cannot be incorporated 
in the Dirac description valid in the long wavelength limit. We can however model the combined effect of these short 
range interactions through suitable boundary conditions. We shall be guided in our search for the suitable boundary 
conditions by the requirement of self-adjointness of Hj. , which is essential for the conservation of probabilities and 
for the unitary time evolution of the system. For gapless graphene, this has already been discussed in 21 j which we 
shall review and elaborate below. 

The operator is symmetric (or Hermitian) [so'l in the domain D{Hr) defined by 

D{Hr) = I *(0) = ^'(O) = 0, absolutely continuous, * e L'^{rdr)}. (10) 

In order to see if the operator Hr is already self-adjoint in D{Hr) we follow the procedure of von Neumann (tI, Isot. 
which requires the analysis of the square-integrable solutions of the equations 

Hl'^±=±i^±, (11) 

where H}. is the adjoint of Hr and are given by 
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Note that H,. and Hi have the same expressions as differential operators although their domains could be different. 
If the Eqs. pT|) have no square-integrable solutions, then the operator Hj. is self-adjoint in the domain D{Hr). 
However, if there are equal number(s) of square integrable solutions of Eqs. pip , then the Dirac operator Hr is 
not self-adjoint for the boundary conditions that define the domain D{Hr) in ITUl) . The latter could happen due to 
the effect of short range interactions induced by the impurity. If that happens, we would have to find a new set of 
boundary conditions, or equivalently a new domain for the radial Dirac operator in which it would be self-adjoint. 
We begin our analysis with Eqs. pT|) . which lead to a set of coupled equations 



du± 
dr 



+ (j/ + i/3 zp 2r)u± - jv± = 0, 



dv± 

r— \- — ip)v± — ju± ~ 

dr 



(13) 



(14) 



for the functions u±, v± appearing in ()12|) . 

In terms of the variable z = —2ikr, equations (fT3l) and ((M]) can be combined to give 



(15) 



Here it is implicitly understood that z — ±2r for k — zLi, respectively. As mentioned before, we have to search for 
square-integrable solutions of Eq. pTj) . We first look at the case with k — +i. The solution of (fT5|) which leads to a 
square integrable solution of (jlip is given by 



v+ = 1/(1^ -iP, 1 + 2v, z) = U{v -if], 1 + 2v, 2r). 
With the help of Ea. p4)) . this result for v+ further leads to 



u+ 

so that we finally have 

F+{r) 
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U{v -iP, 1 + 2v, 2r) - 2rU{v + 1 - 2 + 2v, 2r) 



U{v-il3, l + 2iy, 2r)r''-^ -'^^^^^—^U{iy + l-il3, 2 + 2v, 2r)r''+3 



(16) 



(17) 



(18) 



3 3 
for the radial part of the first component of the vector ^'_|_ spanning the k = deficiency subspace. In deriving 
the expression for we have used the relation C/'(a, 6, z) ~ —all (a + 1,6+ 1, z). 

Next we want to find conditions under which and consequently are square integrable. In the limit 

r — > cx), U{a, b, 2r) ^ r^" and — s- 0, meaning that is square integrable at infinity. As r ^ 0, 

M(a, 6, 2r) 1 and 



\F+\rdr 



'' dr + convergent terms. 



(19) 



Therefore, for the range < J/ < i, F+ is a square integrable function. Similar analysis shows that for < < ^ 
the entire radial wave-function is square integrable. Thus, for k = +i we have a single square integrable solution 
of Eq. (fTT|) . The deficiency index n+ of Hr is defined by the number of linearly independent and square integrable 
functions satisfying Eq. (|lip for k = +i. Our analysis shows that for gapless graphene with a charge impurity, n+ = 1. 

Let us now consider the case when k = —i. In this situation we have z = —2ikr — —2r and a possible solution 
of (fTS]) which leads to a square integrable solution of (|lip is given by 



V- = e/U{l + i^ + iP, l + 2iy, - z) = e-^'^Uil + iy + i(3, 1 + 2iy, 2r). 
By using again Eq. (|14p . this result for w_ yields the following expression for w_ 

u- ^ ^"^^ e~'^'^U{l + v+iP, l + 2i/, 2r)-— e"^'^ U{l+v+il3, l + 2v, 2r) + (l+z^+i/3)C/(2+zy+i/3, 2+2i^, 2r) 
3 3 1 



(20) 



(21) 
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If we make use of the recursive relation U{a, b, z) — aU{a + 1, 6, z) — U{a, b — 1, z) 
be simplified to 



0, last expression for w_ can 



i/3 



2r 



e-^Wil + v + i(3, l + 2v, 2r) ~ —e-^^Uil + v + il3, 2 + 2;/, 2r), 



so that we finally have 

F_(r) = e-' 



r ■ 

j 



U{l + v + il3, l + 2u, 2r)- -r^+^Un+iy + iP, 2 + 2v, 2r) 

j 



(22) 



(23) 



for the radial part of the first component of the vector spanning the k ^ —i deficiency subspace. By going 

through the similar procedure as before, we see that as r — > oo, F^{r) — > 0. In the short distance limit we can again 
show that -F- and the corresponding entire radial wave- function is square integrable when < v < i. Thus, for 
k = —i also, we have a single square integrable solution to Eq. pTj) for the range < < ^ of the parameter ly, 
implying that the deficiency index n_ = 1 . 

According to von Neumann's analysis [sBl, this result n+ = n_ = I, implies that Hr itself is not self-adjoint 
on the domain D{Hr), PH]) . but can be made self-adjoint through a suitable choice of boundary conditions. When 
71+ = 7i_ = 1, the allowed boundary conditions are labelled by a single real quantity, called the self-adjoint extension 
parameter. Our analysis thus shows that radial Dirac operator Hr admits a one parameter family of self-adjoint 
extensions when < v < i. When = n_ = 1, von Neumann's analysis further says that domain in which Hr 
is self-adjoint is given by I)z{Hr) = 'D{Hr) ® '1'+ + e^*t where z £ R (mod 2n) is the self-adjoint extension 
parameter [sol. Physically, the domain T)z{Hr) provides the boundary conditions for which the radial Dirac operator 
for graphene is self-adjoint and we see that the boundary conditions are labelled by the parameter z. 

In subsequent analysis it will be of interest to know the short distance behaviour of functions and _F_ . Their 
behaviour at short distances can be deduced with the help of the expression [29 1 



U (a, 6, z) 



sin nb 

Using (j24|) . in the limit r ^ 0, we have 

TT V + j — i(3 



M{a,b,z) 



r(i + a - b)r{b) 



Mil + a~b,2-b,z) 
r(a)r(2-6) 



1 



sin7r(l + 2i/) j T{~iy - iP)T{l + 2iy) 



(24) 



(25) 



-1v 



_sin7r(2 + 2z/) j r(l + t/ - i/3)r(-2i/) sin7r(l + 2i/) j V{v ~ - 2v) 



and 



TT 



iy + j-i(3 



1 



sin7r(l -I- 2iy) 

TT 1 



j T{l-v + il3)V{l + 2v) 

2-2"^ TT 



(26) 



_smT:{2 + 2v) jT{l + v + il3)T{-2v) sin 7r(l -t- 2j/) j T{1 + v + il3)V{l - 2v) 
This concludes our discussion of the allowed boundary conditions for the radial Dirac operator for the gapless graphene. 



2.3. Phase shifts and S matrix for gapless graphene with a charge impurity 



We now turn to problem of finding the scattering state solutions to Eq.(IT]) for the parameter range < < 



This problem amounts to solving the set of coupled equations ([7]) and ([5]) or, equivalently, finding the solution of the 
2nd order Eq.®. The required solution of Eq.® which leads to physical scattering states has the form 29 1 



V = AiM {v -if], 1 + 2v, z) + Aaz^^^'M {-v -if], 1 - 2iy, z) , 



(27) 



where z = —2ikr. 
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By assuming v in the form (|27|) and with the help of Eq.®, u can be calculated as 



ju^ Aiz^M{v-il3, l + 2u, z) + A2z{~2u)z-'^''-^M {-V - il3, l-2u, z) + A2Z-^''+^ -^M {-v - ip, \-2v, z) 



d_ 
dz 



dz 



(28) 



+ Ai{iy-i(3)M{iy-iP, l + 2iy, z) + A2{v - il3)z-'^'' M {~u ^ if}, l~2v, z) . 
After applying the recursive relation 

aM(a, 6, z) + zM'{a, b, z) = aAf (a + 1, 6, z), (29) 
where prime denotes a derivation with respect to z, we get 

ju ^Ai{iy^ i(3)M (1 + 1/ - 1 + 2iy, z) + A2Z-'^''{-v - if3)M {l-iy-if3, 1 - 2iy, z) . (30) 

This leads to 

F(r) =r''^ 5 Ai'^-^M{l + iy-il3, l + 2v, z) + A2Z-^''^^^^^M [l ~ v - i(3, 1-2;/, z) 
V 3 3 



+ AiM {u - iP, l + 2v, z) + A2Z-^''M{~v-il3, 1 - 2iy, z) 



(31) 



for the radial part of the first component of the wave- function in ([3]) . 

In order that the Hamiltonian provides a unitary evolution, the physical solution obtained above must belong to 
the domain of sclf-adjointness given by 'Dz{Hr). We will enforce this by matching the behaviour of the physical 
wavefunction with that of a typical element of 'Dz{Hr) in the limit r ^ 0. In order to do this, note that as r — > 



F{r) ^AAl + ^) r^-'^ + A2i-2zk)-'- ( 1 + Z^^) ^ 



In the same limit a typical element of the domain 'D^_{Hr) is given by 

^f(r, (j)) = \ (e*4«'+ + e-*^ , 



(32) 



(33) 



where A is a constant. Comparing the upper component of the physical solution to that of an element of Vz{Hr) in 
the limit r — )■ , we get 



Ai = A- 



6 2 



sin7r(l + 2v) yTi-u - i(3)T{l + 2iy) r{l - ly + i/3)r(l + 2iy) J ' 



(34) 



A2{-2ik)-^''{j-iy-if3) = -X- 



sin7r(l + 2i>) 



'e'2 



T{l + v -iP)T{~2u) 



'e'2 



T{v -iP)T{l~2v) 



T{l + v + iP)T{-2v) 



+ + iP) 



V[\ + v + iP)V{\-2v) 



(35) 



We now proceed to find the phase shifts characterizing the scattering process in gapless graphene. These quantities 
are of interest to us because it is possible to extract important informations out of them, which in turn are related to 
various physical properties of graphene such as the electrical conductivity and transport properties. In order to find 
the phase shifts and the scattering matrix, we have to investigate the asymptotic behaviour of the wave- function ^P, 
Eq.Q, or equivalently, of the function F(r) given in the relation (pij) . To do this, we note that the behaviour of the 
confluent hypergeometric function M{a, b, z) in the asymptotic region Re(z) = and Im(z) — > -f oo is given by [2^ 



M{a,b,z) 



m 

r(«) 



l + Oi\z\-') 



m 

T{b-a) 



l + 0{\z\'') 



(36) 
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Due to the fact that we are deahng with the problem where Re(z) = 0, both leading terms in the asymptotic 
expansion (|36p of AI approximately have the contribution of the same order, so that both of them have to be taken 
into account. Knowing that, the behaviour of F{r) when r ^ oo is now given by 



F{r) ^ {-2ik)-''{^i)-'P ( ' + A 



i(fer+/31n2fer) 



The scattering matrix and the corresponding phase shift can be written as 



S{k) EE e^'^(fe)e-^-^^ r 'Tr Z^ ^ 7 r( 2 ^ ^ (38) 

Ai J r(i+i/-i/3) + j r(i-iy-i/3) 

where the ratio A1/A2 can be obtained from Egs. dMl) and (IHS)) and is given by 

4i ^ ~ ^'^j ^^(^ ^ r(CT gT_ (39) 

^2 r(l + 2l/ e't , e-'t ^ ^ 

r(i/-i/3) ^ r(i+i/+»/3) 

It should be emphasized that the above analysis is relevant only when the system parameters are such that < < ^ , 
where v = \J p — fP-. For this range of the system parameters, the radial Dirac operator Hr in Eq. ^ is not self- 
adjoint with the boundary conditions describing the domain D(Hr) (|10p . The allowed set of more general boundary 
conditions for which the radial Dirac operator Hr is self-adjoint is labelled by a real parameter z defined modulo 2tt. 
The above analysis clearly shows that the spectral data including the phase shifts and the scattering matrix explicitly 
depend on the boundary condition through the self-adjoint extension parameter z. For each value of z, we have an 
inequivalent quantum description of the gapless graphene system. This is not surprising as in our framework, the 
boundary conditions (or equivalently the self-adjoint extension parameter z ) encode the effects of the short range 
interactions due to the impurity. Different values of z correspond to different combined effects of the short range 
interactions, leading to inequivalent physical results. It should be noted that the value of the self-adjoint parameter 
z cannot be fixed by the above analysis alone. For any given system of gapless graphene with a charge impurity, the 
value of the parameter z would have to be fixed empirically by measuring physical quantities which depend on the 



scattering data, such as resistivity 28|. We shall make further comments about this in the conclusions. 

We end this Section by noting that the S matrix in (j38p docs not have any poles in the positive imaginary value of 
the energy. This implies that there are no bound states for a gapless graphene in the presence of a subcritical strength 
charge impurity, in agreement with the Klein paradox. 

3. SCATTERING IN GAPPED GRAPHENE WITH A CHARGE IMPURITY 

In this Section we shall discuss the scattering states, phase shifts and the S matrix for a gapped graphene in the 
presence of a charge impurity. We start with a brief review of the gapped graphene system where a charge impurity is 
present. This is followed by a review of the allowed boundary conditions for which the gapped graphene Hamiltonian 



admits a unitary evolution [12| |. Finally we find the scattering data for the gapped graphene system which under 
certain conditions depend on an additional parameter that labels the boundary conditions. We also comment on how 
to determine this parameter empirically. 

3.1. Gapped graphene with a charge impurity 



In this case, the 2D Dirac equation can be written as 



(40) 
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where the Hamiltonian is given by 



a 



H = -i{(Jidx + cr2dy) + ma'i (41) 



r 



and the eigenfunctions are two-component wavefunctions of the form 



(42) 



In the above equations r and <j) are the radial variable and the angle in the x — y plane, respectively. Consider 
the ansatz 

i;i{r) = V^^{TE e-^p''-i{P{p) + Q{p)), (43) 



i;2{r)^^M^ e-^p''-^{P{p)-Q{p)), (44) 



where p — 2'-fr = 2\/m^ — E'^ r and v = ^/ p — a^, j being a half integer. After inserting the ansatz (j43|) . ((44|) 
into we obtain the set of coupled equations 

where Hp defined above denotes the radial Dirac operator. These equations can be decoupled to give 

(fP „ ,dP f aE 

dp'' dp \ 7 



+ (1 + 2^--°)— - (^- — )^ = 0' (46) 



p^ + (l + 2^-p)^- (l + i.- — )g = 0. (47) 



These equations can be solved in terms of confluent hypergeometric functions [29j. We now proceed to discuss the 
allowed boundary conditions for this problem. 

3.2. Boundary conditions for gapped graphene with a charge impurity 

The analysis for the allowed boundary conditions for gapped graphene with a charge impurity was already discussed 



m our previous paper 



12( 1 . Here we review that work for completeness and also for obtaining the formulae that would 
be useful in the subsequent analysis. In order to find the deficiency indices n± for Hp and subsequently for iJ, we 
need to solve equations 

i7^*±=±i5'±, (48) 

where iJ^ has the same differential expression as for H in Eq. (j4ip . although their domains could be different. In 
Ea. (|48p are two-component spinors of the form 



where 



i^i± = \A^e-^p'^-^(F±+g±)(p), (50) 
V;2± = V^^te-^p''-HP±~Q±)ip)- (51) 
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The functions ipi± and ip2± have to be square integrable in i?+ with a measure pdp. To get any further information 
one has to solve for P± and Q± from the equations 

dP+ , ia ^ ( ma \ „ , , 

dp 7± V7± / 

dp 7± V 7± / 



where 7± = \Jrn? + 1 and p = 27±r . These equations can also be decoupled to give 

-—^ + (l + 2iy-p)-^- U 

dp^ dp \ 7± 

(PQ± r, •,dQ± ( ia 



P^ + {l + 2.-pr-^-{.--]P^ = (54) 



+ {l + 2iy- p)^- + u g± = 0, (55) 

dp \ 7±/ 

Let us first look at the deficiency subspace determined by the positive sign in the above expression. Then the required 
solution of Eq. is 

P+-C/(^i.-^,l + 2z.,p^, (56) 

where U denotes a confluent hypergeometric function defined in Eq. (|24p . Since is not independent of P+, it can 
be calculated with the help of the relation (|52p . 



3 ]Q+ = p-rU{y A + 2v,p] + [y ]u[v ,f + 2i.,p . (57) 

7+ / dp \ 1+ J \ 1+J \ 1+ J 

By using the differential recursive relation zU'{a, b, z) + aU (a, 6, z) — a(l + a — b)U (a + 1, 6, z), with prime denoting 
a derivation with respect to z, we get 

'ma\ ^ ( iCt\ f iCt \ rr f . , r, \ /^r,N 

^ --)(-"- -J ^'^''--''-^''''T '''' 

In the limit p — > the functions ((56)) behave as 

P+ ^ a{A+ - B+p-^n, (59) 

g+ — > a(C+ - (60) 

^here a = and 

^+ = r(-;.-^)r(i + 2^.) ^+ " r(i. - i^)r(i - 2^.) ^^^^ 

+ (j-^) r(i - - i^)r(i + 2;.) + (j-^) r(i + 1. - ^)r(i - 2;.) ^ ^ 

are constants depending on the system parameters. We can now easily show that and '02+ are square integrable 



everywhere provided < v < ^ [1^. Thus n+ = 1 for the parameter range < v < ^ 



In a similar way, by analyzing the deficiency subspace characterized by the negative sign in (jiSj) . we obtain 



la 



P- =U + —,l + 2:y,p] , (63) 



ma \ „ la \ la \ la . , , 

J ]Q- = [v+ — ][-v+—]u[l^v+—,l + 2v,p]. (64) 

1+J \ 1+J \ 1+J\ 1+ 
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In addition, in the limit p — > the functions P_ and behave as 



Q- 



where a 

A- 



\tt(1+2v) 



and 



1 



C- = 



1+ ' 



1 



(j - ^) 



From this we see that relations 



^)r(i-2i/) 



7+ 



r(i 



)r(i - 21.) 



(65) 
(66) 



(67) 
(68) 

(69) 



hold. Similar analysis as before shows that n_ = 1 for the parameter range < // < ^ as well. Thus for the gapped 
graphene with a charge impurity, 71+ = n_ = 1 when < i/ < i. Therefore, this system admits a one parameter 
family of self-adjoint extensions for < < i 12| . 



3.3. Phase shifts and 5* matrix for gapped graphene with a charge impurity 

Scattering states correspond to solutions of Eq. lHDl) when E > m. In this case the parameter 7 = %/ — E"^ 
becomes purely imaginary, i.e. 7 = iq, where the real parameter q is defined as g = \/ E^ — m^. Consequently, the 
variable p also becomes purely imaginary, p = 2iqr. 

In order to analyze the r — > cx3 limit of the scattering states, we first express the solutions written in terms of 
the hypergeometric function U in terms of the hypergeometric function AI using Eq. (|24p and then carry out the 
asymptotic expansion p6[) . As was the case for the gapless graphene, here again both leading terms in the asymptotic 
expansion (j36p of M approximately have the contribution of the same order, so that both of them have to be taken 
into account. 

The solutions of Eas. (|15)) . (|T7)) which lead to physical scattering states are found to have the following form 

P{p) = AM + 2v, + Bp-^-'M (^-1. - 1 - 2^, p^ , (70) 

Q(p) =CM(^l + i.-^,l + 2i/,p^ + Dp-^'' M (^1 - - ^,1 - 21., ■ (71) 
In this case the solution for the function ipi, appearing in would look like 

A{q)M (^iy~^^l + 2i.,p^ + B{q)p-^''M (^v~'^,l~2v,p 



^pi{r) ^ Vm + E e ^ p" 



( OiE \ ( aE 

+ C{q)M h — -,l + 2l/,pj +D{q)p-'^''M fl - 1/ - — , 1 - 2j/, p 



(72) 



where we have assumed that the coefficients A{q) and B{q) depend on the real parameter q. However, at this stage, it 
is important to stress that functions P and Q are not independent of each other, but rather they are related through 
the set of coupled equations 

dP , aE . ^ ( ma \ ^ „ ,„„x 

p-r^^"" )^+ J Q = o> (73) 

dp 7 V 7 / 
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do , aE ^ ^ ( ma\ „ 

dp 7 V 7 / 

As a consequence, the constants A{q), B{q), C{q) and D{q) are also not all independent. 
By assuming that P has the form ([TO)) , we can obtain Q as 



(74) 



J 



7 



7 



(75) 



AM' - 1 + 2iy, + Bp-^^M' (^-f ~^,l~2u,p^ + B{-2v)p-^''-^ M (^-v -^,l-2v,p 

where prime denotes the derivation with respect to p. By applying the recursive relation 

aM {a, b, z) + zM'{a, b, z) = aM{a + 1, 5, z) 
(prime denotes derivation with respect to z) we obtain 



(76) 



J - — j g = A{q) iu- — \Mil + iy-—,l + 2:,,p]+ B{q)p- 



2u 



aE\ , /, aE 
-V \M\\-v , 1 - 2i/, p 

7 / V 7 

(77) 



Since Q from ([77]) has to coincide with Q from ([7T|) . (at least up to a constant), we must have 



(78) 



D{q)^B{q)- 



(79) 



As r — > oo, using ([36]) and noting that 7 = ig, the ^1 component of the wave function has the form 



■ipi{r) > Vm + E {2iq) ^ ' 1 



+ Vm + E {2iq) 2+* 1 



Aiq) 



C{q) 



T{l + 2v) 
T{l + 2v) 



r(i - 2v) 

\ q } 



r(l - 2iy) 



(80) 



from which one can identify incoming and outgoing waves and the corresponding amplitudes. The scattering matrix 
and the phase shift can be obtained from the limiting form (jSOp of the wave function as a ratio of its outgoing and 
incoming amplitude. 



r(i-2t/) 



S{q) = e 



{2iqr- 



' C \ * Q J 



(81) 



So far the constants A{q)^ B{q), C{q) and D{q) are not completely specified. All we know about them is that they 
are related through the relations ([78]) . (f79|) and that they depend on the transfer momentum q. To find a further 
relationship between them, it is of interest to look at the short distance limit, r ^ 0, of the function ((72|) 



^^(^) ^ V^^rT:^ ((A(g) + C{q))p''-^ + {B{q) + D{q))p—^) 



(82) 



This expression gives us the function ■01 to lowest order in r. 

We already know that the radial Dirac operator Hp and consequently Hamiltonian iJ, (|¥T|) . admits a one parameter 
family of self-adjoint extension when Q < v < ^. In this case, the domain of self-adjointness of the Hamiltonian ((4T|) 



12 



is given by T>z{H) = 'D{H) © v]/^ 4- e *2 \]/_}. In the limit r — > 0, a typical element of the domain T>z{H) can be 
written as 

*(r,(/.)=c(e***+ + e-*4*_), (83) 



where c is some constant and are square integrable solutions of Eqs. 

After we make use of the relation ([83| and the fact that the coefficients of appropriate powers of r at both sides in 
([83)1 must match, we get the following two conditions 

(27)"^"^ {A{q) + C{q)) Vm + E = ca (v^^TTIe't [A+ + C+) + xA^T^e-^t (A+ + (7+)) (84) 

(27)-''"5 {B{q) + D{q)) Vm + E = -~ca (^Jn^ie'^i {B+ + D+) + V^^^ie-'^ {B+ + (27+)-''-^ (85) 

where z is the self-adjoint extension parameter and all other quantities appearing in the above relations are defined 
in the previous section. The last two equations yield 



A{q) + C{q) _ Vm + ^e't {A+ + C+) + - ^e-'2 {A+ + C+) 2u ,^ .-2u _ acos(gi + f) .^2. 
B{q)+D{q) V^^ie'^(B+ + D+) + V^^ie-'^{B+ + D+y ^ ^' 6 cos(02 + §) ^ ^ ^' ' 

(86) 



where we have defined Vm + i{A+ + C+) = ^ie**i and Vm + i{B+ +£»+) = ^je'^". Using ([75]), (HH) and the 
S matrix in (jSip can be written as 

Cicos(ei + f) ^21. N -21.1+^ r(i+2t/) , r(i-2i.) 
2,^ -g.cos(e.+f)(27+) (27) T+^^ir(i+.+.^) +^2r(i_,+,^) 

where 

-1 = ^-—^' "^=s(^-7^- ^''^ 

The expression in (|87p gives the 5* matrix for gapped graphene for the parameter range < ly < ^ . For this range 
of V, the appropriate boundary conditions for which the Hamiltonian (|4ip is self-adjoint and the corresponding time 
evolution unitary requires the introduction of an additional real self-adjoint extension parameter z, which labels the 
allowed boundary conditions. The phase shifts and the S matrix depend explicitly on z. For each value of z (mod 
27r) , we have an inequivalent set of the scattering data. The above analysis by itself cannot determine which value 
of z would be realized in a given physical situation, which must be determined empirically. We shall make further 
comments about this in the conclusions. 

We end this Section by noting that we can recover the bound states found in 12] from the poles of the S matrix 
(|87p in the positive imaginary axis. In the bound state sector m > E, hy setting q — ip, p — y/m^ — E^, as some 
arbitrary pole of the scattering matrix (|87p . we obtain the condition determining the bound state spectrum as 



^2^^r(i-2.)r(i + .-^^) + ^ + acos( 



ljTil + 2u)T{l-u-^^^)\j-^+J^^I^j 6cos( 



(89) 



This equation determines the bound state energies in gapped graphene in the presence of a charge impurity, for the 
parameter range < < ^ [12 1. 



4. CONCLUSIONS 



In this paper we have analyzed the scattering sectors of the Hamiltonians for the gapless and gapped graphene in the 
presence of a charge impurity. We have shown that for a certain range of the system parameters, the corresponding 
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Hamiltonians are not self-adjoint and therefore do not generate unitary time evolution with the usual boundary 
conditions. However, they can be made self-adjoint by a suitable choice of the boundary conditions, which requires 
the introduction of a real parameter in the problem. This self-adjoint extension parameter z is a real number (defined 
modulo 27r ) and for each value of 2;, we have an inequivalent quantization and spectral data for the system. In 
particular, the scattering state wave functions, the phase shifts and the S matrices explicitly depend on the parameter 
z. This parameter however cannot be determined by theory alone, and has to be fixed empirically. 

There are a variety of possibilities to fix this parameter empirically. First, the scattering matrix and the differential 
scattering cross-section determines the resistivity of the system [28|. A measurement of the resistivity would be a 
possible way to fix the parameter z. Second, it is well known that Friedel sum rule is affected by the self-adjoint 
extension [31|, as it depends on the scattering phase shifts which are functions of the self-adjoint extension parameter. 
Thus the Friedel sum rule for graphene [32| would also depend on the parameter z. It follows that the integrated 
density of states would be z dependent and the precise nature of this dependence which can be used to fix the value 
of z is under investigation. 

Our analysis of the scattering sector also yields information about the bound states by looking at the poles of the 
S matrix. Here we have recovered the results for the bound state spectrum which has already appeared elsewhere 
0, 0, [13] ■ Measurement of local density of states (LDOS) using scanning tunneling spectroscopy can also yield 
information about the value of z in a given graphene system. It should be noted that beyond the parameter range for 
which the self-adjoint extension is required, the established analyses of the scattering sector is unchanged. 

The analysis presented here can be extended to a variety of physical effects in graphene. These include the analysis 
of the impact of generalized boundary conditions on screening effect in graphene 33|, |3J, |35|, l36|i l32| and on the study 
of impurities in bilayer graphene [i^l ■ These are presently under investigation. 
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